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Abstract 

Results of numerical simulations obtained by a staggered finite differ¬ 
ence scheme together with an efficient immersed boundary method are pre¬ 
sented to understand the effects of the shape of three-dimensional obstacles 
on the transition of a boundary layer from a laminar to a turbulent regime. 
Fully resolved Direct Numerical Simulations (DNS), highlight that the closer 
to the obstacle the symmetry is disrupted the smaller is the transitional 
Reynolds number. It has been also found that the transition can not be re¬ 
lated to the critical R *. = U^k/v used in the past. The simulations highlight 
the differences between wake and inflectional instabilities, proving that two- 
dimensional tripping devices are more efficient in promoting the transition. 
Simulations at high Reynolds number demonstrate that the reproduction of 
a real experiment with a solid obstacle at the inlet is an efficient tool to gen¬ 
erate numerical data bases for understanding the physics of boundary layers. 
The quality of the numerical method to fully resolve the small scales, that 
is one ingredient for a DNS was shown by a comparison of the exponential 
range of the velocity spectra, in Kolmogorov units, with those for isotropic 
turbulence. The good comparison reinforces the idea of local isotropy at the 
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smallest scales. 
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1. Introduction 


The laminar-turbulent transition of boundary layers is of large impor¬ 
tance in several practical applications and in particular for the aeronautical 
industry. Naturally it occurs when external disturbances affect the near wall 
region, generating turbulent spots (Emmons [(J) occurring randomly and 
growing uniformly independently of one another as they are swept down¬ 
stream by the flow. More efficiently, the transition can be promoted by in¬ 
serting solid obstacles on a smooth wall at a certain distance from the leading 
edge. Therefore the shape of the elements may play an effect on the location 
at which the flow becomes turbulent. The tripping device can be obtained 
by a series of obstacles widely separated (thee-dimensional) or by a wire in 
the whole span (two-dimensional). Different scenarios can be achieved, for 
the former configuration the transition can be classified as a wake instability, 
which is associated to the generation of vorticity layers at the obstacle walls 
which evolve in the streamwise direction. For disturbances in the whole span 
of the plate an inflectional point in the velocity profile inside the recircu¬ 
lating region behind the obstacle appear. Therefore the small disturbances 
grow rapidly and a turbulent flow is generated in the whole boundary layer. 
Klebanoff et al. 0 wrote an exhaustive paper on the transition of boundary 
layers induced by solid elements, claiming that two-dimensional obstacles are 
more efficient to excite three-dimensional disturbances than the oscillating 
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ribbon previously used (Klebanoff et al. 16j). The Klebanoff study had 


such impact that this type of laminar-turbulent transition was defined as 
Klebanoff transition. 

Klebanoff et al. [ijJ used hemispherical obstacles with an aspect ratio 
equal to 1 and height k smaller than the boundary layer thickness 5. They 
compared the values of the transition Reynolds number (R/. = Ukk/u, where 
Uk is the velocity of the undisturbed flow at the top of the element) with 
those obtained in other experimental setup. Large differences were found, 
for instance their value equal to 325 was different than those in Tani et 


al. 


32], varying between 600 and 1000. Dryden [5] emphasized that, once 


the transition occurred, a further increase of the Reynolds number had the 
effect of moving the location of the transition upstream. The reasons of the 
differences in the critical R.p. depend on the nature of the laminar-turbulent 
transition. For instance for isolated obstacles it occurs at a certain distance 
from the obstacle rather smoothly, and each scholar used different criteria to 
establish when the turbulent flow is assessed. In the discussion of the results 
Klebanoff et al. jl^] claimed that the shape of the obstacle plays a role and 
that the data in literature were limited. For us this assertion is still true and 
hence it has been decided to give a contribution to understand the effects of 
the shape of the obstacles by numerical simulations. 

The Klebanoff’s transition was reconsidered experimentally and numer¬ 


ically by Bake et al. j^j]. However, the DNS was performed by taking the 
disturbances from the experiments and it was observed that the vortical 
structures responsible for the transition consisted in series of A-shaped vor¬ 
tices. Ergin & White jii], analysing the experimental results, related the 
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transitional mechanism to the shape of the optimal disturbances produc¬ 
ing the largest growth over a specified distance or time. They claimed that 
these disturbances are similar to those produced by three-dimensional ob¬ 


stacles. Flow visualizations as that by Acarlar & Smith [1] have shown that 
three-dimensional obstacles in a laminar boundary layer generate horseshoe 
vortices wrapping around the upstream side of the obstacle with the two legs 
downstream. Flow visualizations are qualitative and quantitative measure¬ 
ments are difficult, in particular to have insights on the effect of the shape of 
the obstacle. Acarlar & Smith [1] for clear flow visualizations used obstacles 
with k/S = 1 different from those typical of the tripping devices and differ¬ 
ent from those used by Klebanoff et al. [17J]. Therefore, they got different 
profiles for the streamwise velocity fluctuations. In discussing the reason of 


the differences and on the possibility to find a scaling of the transition, K 


banoff et al. [17|] were asking for further investigations. Ergin & White 


al 


designed their experiments to understand the flow behavior in the Kendall’s 


near and far wake (Kendall [15|). In the near wake region the energy of the 
disturbance decreases; in the far wake, depending on the Reynolds number, 
the growth is similar to the optimal growth predicted by the theory. Ergin 
& White [9] used an array of cylindrical elements separated by a distance 
such that the elements can be considered isolated. They did not change the 
shape of the obstacles and they still used the critical Rk as a parameter to 
establish when the flow becomes turbulent. 

Our opinion is that Rk can not mark the transition from laminar to tur¬ 
bulent flow. The Rk should be replaced by a quantity related to the velocity 
disturbance generated by the solid obstacle. In wall bounded turbulent flows 


4 






the normal to the wall velocity fluctuations (u 2 ) are those producing the 
turbulent Reynolds stress (u\u 2 ) characteristic of the turbulent regime (the 
index i — 1 indicates streamwise, i = 2 normal and i = 3 spanwise di¬ 
rections, angular brackets () averages in the homogeneous directions and in 
time). These fluctuations produce also (u 2 ), the only normal stress entering 
in the mean momentum equations, for channel and pipes. For boundary lay¬ 
ers other stresses appear, but the relative terms are one order of magnitude 
smaller. In wall bounded flows the importance of ( u 2 ) was not emphasized 
in the early experiments, designed to understand the physics of the near wall 
turbulence, for the difficulties to measure U 2 , close to a wall, by hot wires. 
In numerical simulations all the quantities are available and in addition the 
boundary conditions at the surface can be arbitrarily varied. For instance, 
in turbulent channels, Orlandi et al. 23] demonstrated that, the normal ve¬ 
locity distribution on the plane of the crests is the driving mechanism for the 
modifications of the near wall structures. The preliminary results suggested 
that, a new parametrisation for rough flows could be obtained by u' 2 \k, with 
u[ = (uf) 1 / 2 (|fc ndicates values at the plane of the crests). A continuous 
transition between smooth (u' 2 \k = 0) and rough walls (u 2 \k ^ 0) was ob¬ 
served. The importance of this quantity has been further emphasized by 


Orlandi 


20] performing DNS of transitional channels with different type of 


solid disturbances in one wall of the channel. It was observed that the flow in 
the channel becomes turbulent when u 2 1£ > 0.2 (the superscript + indicates 
wall units). However when this threshold is achieved, the smallest value for 
turbulent flows is u 2 \= 0.6. By increasing the Reynolds number there is a 
growth of this stress tending to a value approximately equal to 1, which is 
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also the value of the maximum v! 2 in smooth pipes or channels. 


Bake et al. [2| introduced normal velocity disturbances at the wall derived 
by laboratory experiments designed ad hoc. From the previous arguments 
a similarity between a solid obstacle and a distribution of w 2 at the wall, 
for instance a round jet, does exist. Therefore for rough boundary layers u' 2 
could be a quantity whose value when overcomes a threshold value detect the 
laminar-turbulence transition. Regarding the transition this quantity can not 
replace the R k criterium because it can not be apriori extimated, but quali- 
tatitevely it can help to design the optimal shape of the obstacle producing 
the highest value at the top, this is out of the scope of the present paper. 
The strategy to change the shape of the obstacle helps to demonstrate that 
R k can not mark the transition. In fact for certain geometries the transition 
may occur and for others no, even if R k does not largely varies. A further 
interest is to investigate whether the numerical simulation of a laminar inlet 
flow past a solid obstacle can be an affordable procedure to get turbulent 
boundary layers at high Reynolds number. To reach these goals it is nec¬ 
essary to understand which tripping device gives the fastest growth and to 
have a quantitative measure of the transitional distance necessary to get a 
fully turbulent regime. The first DNS of boundary layers by Spalart 31] used 


a pseudo-spectral numerical method together with a fringe method to recy¬ 
cle the output flow to the inlet. More sophisticated recycling methods were 
developed as reported by Pirozzoli & Bernardini [25] in DNS of compressible 


turbulent boundary layers. Simens et al. 


30] asserted that with a distance of 


300 momentum initial thickness the influence of the artificial inlet conditions 
are eliminated. Lee et al. 18] used similar inlet conditions in presence of 
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rough walls. In these conditions the disturbances, continuously emanating 
from the interior of the rough surfaces, reduce the distance to achieve a fully 


turbulent rough boundary layer. Wu & Moin 


33] simulated the whole tran¬ 


sitional flow by adding to the potential flow patches of isotropic turbulence. 
This disturbance can be considered an attempt to reproduce the conditions 
in wind tunnels similar to those by Emmons 6], where turbulence grows as 


random organised turbulent spots. More recently Schlatter et al. 


271 ] used 


a low-amplitude trip force acting in the wall-normal direction. In this way 
they had a rapid laminar-turbulent transition close to the inlet. The au¬ 
thors did not furnish the velocity distributions so it is difficult to connect the 
disturbance to a real tripping device. 

With the exception of the DNS of rough boundary layers (Lee et al. [l3 | ) 
the numerical methods in the mentioned simulations did not allow to repro¬ 


duce the flow interaction with so 
this paper (Orlandi & Leonardi 


id obstacles. The numerical method used in 


21]) is based on a second order finite differ¬ 


ence method and on an immersed boundary technique (IBM) to reproduce 


bodies of any shape. The method was 
ticular it was shown by Orlandi et al. 


ested in previous papers and in par- 


22] that the pressure distribution on 


the rod-shaped elements was in a very good agreement with those measured 
by Furuya et al. 11]. In addition a combined experimental-numerical study 
(Burattiui et al. ^41) was performed to validate the DNS, with square bars 
with w/k = 3 (w is the spacing between two elements). The basic numerical 
method requires small modifications to deal with boundary layers. Radiative 
outlet boundary conditions allow to the vortical structures to exit smoothly 
from the computational domain. The modifications in time and space of 
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the three components of the vorticity fields, near the solid element, allows 
to understand why at the same Reynolds number for some geometries are 
amplified and for others are not. A comparison between coarse and refined 
simulations (not presented) suggests that the inaccurate representation of the 
surface acts as disturbances reducing the critical Reynolds number, therefore 
small inperfections in a solid obstacle while do not change Rk would enhance 
the transition, this is a further reason of the inadequacy of Rk to mark the 
transition. 


2. Numerical method 


The basic numerical method consists of a second order finite difference 
scheme with staggered velocities, which in the inviscid limit conserves total 
energy. The method for the channel is described in Orlandi 19] where the 


global conservation properties in presence of non-uniform grids in the nor¬ 
mal directions are presented. The boundary layers simulations require inlet, 
outlet and initial conditions. The radiative conditions at the outlet (Pauley 


et al. 


24]) allow the flow vortical structures to exit from the computational 


domain without producing disturbances, which propagating backwards can 
affect the flow. In these circumstances, a cos FFT is necessary to have a 
direct flow solver for the elliptic equations necessary to have a solenoidal 
velocity field. At the inlet a Blasius velocity profile is imposed by fixing the 
Reynolds number and the distance from the leading edge. The theory of 
laminar boundary layers allows to derive at t — 0 the streamwise and normal 
velocity profiles in the entire computational domain. This two-dimensional 
flow is replicated in the whole spanwise direction, and small random distur- 




bances are added at the inlet to trigger the spanwise disturbances, which, 
interacting with the large scale created by the obstacles at the supercritical 
Reynolds number, produce a turbulent flow. 

The Navier-Stokes and the continuity equation in dimensionless form, are 


du 


duiUj 

dt dxn 


dp 1 d 2 Ui dui 

dxi Re dx 2 ’ dxi 


= 0 . 


( 1 ) 


where p is the dimensionless pressure. The Reynolds number is defined as 
Re = U e S/u , being U e the constant free-stream velocity and S the boundary 
layer thickness at the inlet section (usually indicated as 5 99 ) . 

The solid element is located at a distance Xq from the inlet, whose ef¬ 
fect on the flow is reproduced by an immersed boundary technique. The 
efficiency of the method together with the MPI (Message Parallel Interface) 
instructions allow to use a sufficient number of grid points to describe the 
:>ody surface. The IBM method differs from that developed by Fadlun et al. 


le surface were evaluated 


lo| . where the velocities at the nearest points to t 
by interpolations. Indeed Orlandi & Leonardi 21], after the evaluation of 


the right hand side of the discretized equations in the regular grid, correct 
the viscous terms, in the explicit and implicit step. For the flows here consid¬ 
ered where the obstacles do not move, the metrics coefficients at the nearest 
points are calculated only once initially. 


3. Results 

3.1. Low and medium Reynolds number 

A first set of simulations were performed on a coarse grid with 512 x 128 
points in the streamwise and spanwise directions respectively. In the normal 
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direction a non-uniform grid was used with 50 grid points for 0 < x 2 < 1 and 
50 for 1 < x 2 < 10. The three-dimensional solid obstacles have aspect ratio 
equal to 1 and a dimensionless height k = 0.25. The computational box, for 
the coarse simulations, has Li = 25.6 and L 3 = 3.2. With these parameters 
the solid elements are described by 5 points in x\ and x 3 and 25 points in x 2 . 
Such coarse grid was chosen on the purpose to investigate the dependence of 
the results on the shape of the element. Three-dimensional solid disturbances 
with cylindrical, cubic and wedge shape were used. The experimental results 
described by Klebanoff et al. [ 17 J showed that two-dimensional obstacles 
behave differently than three-dimensional obstacles. Therefore an additional 
set of simulations with a square rod, located at the same distance from the 
inlet (x 0 = 5) were performed. 

These simulations showed that for the cylinder the transition Reynolds 
number was higher than that for the cubic obstacle. The latter did not 
largely differ from that for the wedge, ft was also observed that for square 
bars the critical Reynolds number was much smaller than that for the three- 
dimensional obstacles. From these results it was decided to look at the 
transition mechanism for the cubic obstacle by varying the Re number in a 
well resolved simulations with 1024 x 160 x 192 grid points in a domain with 
the same size above cited. The obstacles is, then, descibed by 10 x 15 grid 
points in x\ and x 3 . The occurrence of transition was detected by monitoring 
the time evolution of of u> 3 at x 2 = 0.0041 and x 3 = 1.6. 

As it was previously mentioned, an unsteady simulation evolves towards 
a stationary condition, that depending on the value of Re, remains laminar 
or becomes turbulent. In both cases the transient time necessary to the 
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a) 



Figure 1: Streamwise profiles of to 3 versus rj = (x — Xo)/k at y/S = 0.0041 at t = 
10,20,30,40,50 and 60 for a) the cube at Re = 5500 t is indicated in the figures, b) 
the cube at t = 60 and diffcrrent Re 10 3 as in the inset. 
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disturbance, produced by the obstacle, to reach the outlet boundary, depends 
on L\. Figure [Ur shows the trends at Re = 5500, where upstream of the 
obstacle is found a high positive peak, due to the formation, of the opposite 
sign vorticity characterising the head of the horseshoe vortex. Behind and 
close to the obstacle the other positive peak is due to the three-dimensional 
separation bubble, which does not depend on time, and weakly is affected 
by Re. When the flow reattaches the wall vorticity varies with Re and 
with time. The profiles at low Re show that the negative bump, while is 
convected downstream, decreases its amplitude for viscous effects, and the 
disturbance does not reach the outlet section. Figure [Ur at Re = 5500 shows 
that after the short transient, ending at 10 < t < 20, the disturbance travels 
downstream without an appreciable reduction for viscous effects. Figure [U> 
with the profiles at different Re show at Re = 4000 the behavior, similar 
to that at Re = 3000, suggesting a critical R cr close to 4500, which leads 
to Rk = 400 not too different from that observed in the experiments. The 
transition starts near the outlet section and this figure demonstrates that the 
simulations reproduce the Dryden jfj] observations of the upstream movement 
of the transition location with the increase of Re. The zoom of the region 
near the obstacle (figure [Ud) emphasizes that, only for Re > 4000, in front 
of the obstacle is negative indicating the presence of an intense horseshoe 
vortex. 

To have a better understanding of the influence of the shape of the vortical 
structures more shapes have been considered at Re = 5500. In figure [2] the 
contour plots of U\. in a small region surrounding the body, allow to get an 
idea of the accurate reproduction of the flow near the body. In addition these 
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contours depict the shape of the obstacles considered: a) cylinder (CY), 
b) cube (CU), c) a wedge pointing upstream ( WU ), d) a wedge pointing 
downstream (WD), e) half cylinder with straight wall upstream (SU) and f) 
with the straight wall downstream ( SD ). The contours at a distance from the 
wall X 2 = 0.15 highlight the shape of the separation bubble. The head of the 
horseshoe vortex does not reach this height then does not appear in figure [2j 
The square bar has been, also, considered (SQ), but for this two-dimensional 
obstacle the geometry is simple and contour plots of u\ are not given. To 
emphasize the shape of the separated region, the Aui, for the negative values, 
has been reduced ten times. In figure [2] the smooth contours, and, hence, 
the absence of numerical oscillations stress the accuracy of the immersed 
boundary technique to reproduce the effects of solid walls. The different 
amplitude of the perturbations, in the region close and downstream of the 
separation bubble is a first indication that for the SU and SD obstacles there 
is the transition. In addition it can be speculated that for SD the transition 
starts closer to the obstacles. Similar perturbations are observed for WD, 
having a smaller amplitude and being more symmetric the transition moves 
downstream. For the the other three geometries no any idea can be drawn 
by these images. 

To have more quantitative results the spanwise averaged Cf — 2 S w /Re 
(Sy } = j^\ w ) has been evaluated by saving eleven fields from t = 60 to t = 70 
separated by one time unit. It can be argued that a more interesting quantity 
should be the u> 3 \ w at = 1.6 plotted in figure la with streamwise oscilla¬ 
tions in the turbulent region. To have smooth profiles the simulations should 
evolve for a long time. To detect the difference in the growth between the 
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d) e) f) 

Figure 2: Sti'eamwise velocity profiles at y/k = 0.6 for the three-dimensional obstacles 
simulations described in a grid with 10 points in xi 15 in £3 and 40 in X 2 ■ The height of 
the obstacle is k = 0.25; red dashed positive, blue solid negative and Aui = 0.025 for the 
red and Aui = 0.0025 for the blue lines. 
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laminar and the turbulent regime by smoother curves a further averaging in 
Xs may help. For the cylinder in figure [3}i there is no transition. For the cube 
the transition is not close to the obstacle. To understand the reasons of the 
transition the surface contours of u 3 and oj 2 may help, as it is discussed later 
on. In figure [2fc> the flow is symmetric with respect to the centerline, which 
is an indicator of the absence of turbulence formation within a reasonable 
distance from the obstacle. A good level of symmetry is also obtained for the 
WU obstacle, where thin layers of high velocity gradients form (figure O:) • 
These being more unstable, produce the transition at a reasonable distance 
from the obstacles, as it results in figure [3]a. The transition is close to the 
obstacle for the other three shapes and to stress the differences, in figure [3jo 
the enlargement of the region with the Cf growth confirms the discussion 
derived from figure [2j The different growth is appreciated only at Re close 
to R cr . At high Reynolds number these differences disappear, and, later on, 
is shown that for the cube, the fully turbulent Cf tends to that of a two- 
dimensional square bar. The interesting result of figure 3 is that the trend 
of the three-dimensional obstacles is different from that of two-dimensional 
obstacles. The transition for the latter flows is sharp, indicating a different 
kind of instability. The inflectional point of the U\ profile, in the separation 
bubble, causes the growth of the disturbances in the whole span. For the 
three-dimensional obstacles the disturbances are the vortical structures gen¬ 
erated by the side walls of the obstacles, producing other patches of vorticity 
which spread laterally to fill-up the entire region. Later on flow visualizations 
are used to describe the different physics. 

For the three-dimensional obstacles R *. = U^k/u (14 is the value of U\ at 
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a) b) 


Figure 3: Streamwise profiles of Cf at Re = 5500; a) in the whole length, b) in a region 
behind the obstacle to see the location of the transition. 


the center of the obstacle and at X 2 = k) has been evaluated. Flow visual¬ 
izations show small differences of u i among the different type of obstacles, 
with Rk varying between 474 and 512, the same range found by Klebanoff et 
al. jl3] in several experiments. Since for some obstacle the transition occurs 
and for others does not, it follows that Rk is not the appropriate quantity 
to establish the transition. A better quantity then should be found. Orlandi 
[20] in two-dimensional channel demonstrated that u 2 at the edge of the ob¬ 
stacle produces a u' 2 \k varying in time. After a transient u' 2 \k was rreaching 
a steady or unsteady condition depending on the Re number. In the latter 
case, charcteristic of a fully turbulent flow, the averaged u 2 in the two ho¬ 
mogeneous directions was greater than a threshold value equal to 0.2. For 
boundary layers u 2 \k is evaluated by averaging in time and in x 3 . It varies 
in X\ then it may indicate whether and where the transition occurs. Figure 
|4jr shows that in the laminar region, before the obstacle, with the small dis¬ 
turbances added u 2 \k is practically negligible. At the obstacle a large jump, 
related to the thin shear layer forming on the top of the obstacle appears. For 
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a) 


b) 


Figure 4: Streamwise profiles of a) u' 2 \k at Re = 5500; b) u 2 1£ 


viscous effects u' 2 \k decreases downstream and for certain obstacles it grows 
again to lead to the turbulent regime. For the WU obstacle this happens at 
rj ~ 60, and, for CU, it may occur at 77 > 80. Figure [4ja, however, stresses 
the improbability to have a turbulent regime for CY. For the other three- 
dimensional obstacles a growth similar to that for WU is predicted, which 
at 77 60 leads to a, u' 2 \k greater than that for the two-dimensional obstacle. 

Numerical and experimental investigations, by different research groups, 
are currently devoted to understand, in channels, pipes and boundary-layers, 
the reasons of the different trends of the maxima of the normal stresses 
with the increase of Re. The common view is that for channels and pipes 
the dependence is weak and the value is close to 1 , for boundary layers the 
dependence is stronger and the values are slightly higher ( « 1.2) see Jimenez, 
& Hoyas [14J]. Figure 0)3 shows that, at this low Re , despite the different 
trend, for the three- and the two-dimensional obstacles u 2 tend towards the 
value typical of boundary layers. For the other geometries, generating small 
disturbances, u 2 decreases for 77 > 5, then if it increases again and it reaches 
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Figure 5: Surface contours of 1 x 3 = —2.5 yellow and L 03 = 2.5 red superimposed to 
w 2 = —2 cyan and w 2 = 2 magenta a) cube Re = 5500, a) SU at Re = 5500. 


u' 2 \k ~ 0.2 it can be asserted that transition to turbulence occurs. The 


value 0.2 is close to the threshold value found by Orlandi 


20 ] in transitional 


channels. It is important to clarify that for boundary layers over smooth 
walls encountering a single obstacle u' 2 \^ varies in x\ therefore it can not be 
considered a threshold parameter to detect the transition. The result here 
found could help in the design of the shape of the obstacle, which should 
produce the highest value for u' 2 \k- 

Flow visualizations of the vorticity field surrounding the obstacle may 
explain the different behaviors discussed for the three-dimensional obstacles. 
The CU and the SU cases are analysed, the first does not lead to transition 
near the obstacle, the latter does. The visualization in figure [5k and figure 
[5b emphasizes the important role of 0 J 2 , related to the spanwise gradients of 
u\ in presence of three-dimensional obstacles. Figure [5k does not produce 
a sinuous instability on UJ 2 , the thin uj 2 layers roll-up remaining stable for 
a long distance. The head of the horses shoe vortex (red structure) at the 
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Figure 6: Surface contours of surface contour of A 2 = 1 a) cube Re = 5500, a) SU at 
Re = 5500. 


bottom of the front face does not play a role. This figure leads to conclude 
that in order to promote the transition the shape of the obstacle should 
create vortical structures with a strong tendency to become unstable. Figure 
[2] showed that for the half cylinders (SU) the thin regions with a clustering 
of Mi isocontours are curved and consequently unstable. The comparison 
between figure Efc> and figure EJr shows equal vorticity layers in front of the 
obstacle. The circular shape of the side wall promotes the instability and the 
breakdown of the LO 2 and CO 3 layers. Flow structures are, usually, visualised 
by the A 2 surface contours (the imaginary part of the complex conjugate 
eigenvalue of the velocity gradient tensor), in figure Eh and figure Eh this 
quantity highlights the horseshoe vortex and a series of hairpins downstream 
the obstacle for the SU, while for the CU there is a similar horseshoe vortex 
and two circular vortices behind the obstacle. These structures have been 


A 


visualised and studied by Acarlar & Smith |1|. More recently a forest of 
hairpins were visualised in the DNS of transitional boundary layers (Wu & 
331]). The A 2 detects the tubular structures which are the structures 


Moin 


19 




Figure 7: Surface contours for the case SQ at Re = 5500: a) W 3 = —2.5 yellow and 
0 J 3 = 2.5, b) u>i = —2 blue and u> 1 = 2 green c) 0 J 2 = —2 cyan and u >2 = 2 magenta d) 
A 2 = 2. 
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less important to drive the instability and the energy cascade characteristic 
of turbulent flows. Figure [5] together with figure [ 6 ] indicate that the rollup of 
the u >2 layers contributes to the legs, the rollup of the 003 sheets to the head 
of the hairpins. 

These visualizations are useful to highlight the differences among the 
three- and the two-dimensional obstacles. A square bar creates a thin vortic- 
ity layer at its top and a recirculating bubble below it. Figure^ shows the 
formation of the thin (U 3 layer bending and breaking at a certain X\ . For the 
effect of the other two vorticity components the CJ 3 layers form structures of 
circular shape. For two-dimensional obstacles and 0 J 2 are similar (figure 
Uh and figure^). These components are generated as fast growing insta¬ 


bilities due to the inflectional point in the profiles of Ui{x 2 ). The 


hairpins depicted in figure 0:1 is similar to that in Wu & Moin 33]. From 


'orest of 


these visualizations the reasons of the sharp growth of the Cf in figure 4 
are understood. The forest of hairpins is observed also by three-dimensional 
obstacles, but at large distances from the solid element. To investigate the 
analogies and the differences between two- and three-dimensional obstacles 
simulations at a higher Re are presented. 


3.2. High Reynolds number 

Having demonstrated that the shape of the obstacle affects the transition, 
it is worth to investigate whether the simulations with the tripping devices 
can produce turbulent boundary layers at a reasonable high Reynolds num¬ 
ber. In the last decade the possibility to use clusters with a large number 
of processors allowed to get DNS competing with laboratory experiments. 
Hoyas & Jimenez [ljj], for plane channels, had results up to R r = 2000. This 
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Re is smaller than those in laboratory experiments, nevertheless for the dif¬ 
ficulty to have accurate measurements at R T = 2000 the results were useful 
to the fluid dynamics community. In particular for the scholars interested to 
the u' 2 i to the vorticity rms, and to high order statistics. Regards boundary 
layers, the largest R T was simulated by Schlatter et al. [2t| at Rg = 4300 (Rg 
is the Reynolds number based on the momentum thickness) and results at 
higher Rg , in progress, were presented at conferences. For boundary layers it 
is easier to reach high Rg in laboratories, as reported in the review by Srnits 


et al. 


29] . However, the numerical simulations helped to understand the 


reasons of certain erroneous conclusions due to the limitations of the size of 
the probes. The present simulations reached Rg ~ 1500, and, to our knowl¬ 
edge are the first one to reproduce the turbulent boundary layers generated 
by tripping devices. 

Simulations for a square bar and for a cube were performed at Re = 7500 
in a domain with Li = 51.2 and Ni = 1024; L 3 was the same as in the 
previous section but the number of points in Xs were doubled. In the normal 
direction the resolution was the same, in fact with 160 points, close to the 
outlet section, the first grid point is at y + = 0.845. The comparison between 
the square bar and the cube showed that, even for the square bar, only for 
certain statistics, a fully turbulent condition was achieved. For the square bar 
a further simulation more resolved in x\ and in a longer domain (L\ = 68.2) 
(Ni = 2048) was performed. In the normal directions the better resolution 
was obtained with L 2 = 9 and 192 points. In wall units the resolution 
in the streamwise and spanwise directions, in the fully turbulent region, is 
respectively equal to 10.5 and 2.7. 


22 




Figure 8: a) Shape factor H\ -j, b) friction factor Cf as function of Rg , the symbols are 


the data in Schlatter & Orlu 
insert. 


281 ] and the lines the present simulations as indicated in the 


Schlatter & Orlu |28| compared the variations of the shape factor H \2 
and Cf with Rg with those by other simulations, showing discrepancies but 
generally the same trend. They reported the values of the fully turbulent 
regime. The present simulations, including the transition, highlight the trend 
towards the fully turbulent regime. Figure Eb reproduces the decrease of H 12 
from the laminar to the turbulent value, which reaches the values of Schlatter 


& Orlu 


at different Rg. For the three-dimensional obstacles the values 


with the cube at Re = 7500 show that the trend for the cube differs from that, 
at the same Re, of the square bar. For the latter geometry small differences 
in the transitional regime are observed, but all the simulations in the fully 
turbulent regime give H 12 in good agreement with that of Schlatter & Orlu 
(28 1 . Figure Eb shows a similar agreement for the Cf. Here the negative 
values for Rg account for the separation bubble. For the square bars at 
Re = 7500 small differences in the Cf are encountered in the transitional 
regime, but independently from the resolution (u/ in the insight of figure Eb 
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indicates the more resolved) in the fully turbulent regime the agreement with 


Schlatter & Orlu 


data extends for a wide range of Re values. 


As previously mentioned, the reason of the different trends to the turbu¬ 
lent regime for two- and three-dimensional obstacles is linked to the different 
type of instability. This is qualitatively described by the velocity contours at 
different distances from the obstacle. Figure [9^ shows that the no-slip condi¬ 
tion on the side of the cube produces two low-speed streaks, better depicted 
by the u >2 contours in a horizontal plane at a distance x 2 = 0.05, in wall units 
approximately equal to 17 (figure 9e). Further downstream other vortex cou¬ 
ples form, which imply an array of low and high speed streaks. In the first 
two x 2 — x 3 planes, in figure 9a, the undulation at the edge of the boundary 
layer is almost negligible, indicating small fluctuations of u 2 , quantified by 
U 2 contours in figure 9c. These disturbances at the fourth and fifth sections 
are enough strong to deform the U\ contours at the edge of the boundary 
layer in figure 9a. The cn 2 contours in figure 9e show the end of the spreading 
at the fifth section and the moderate U 2 ejections, until there, are depicted 
in figure 9c. This figure, further downstream, shows coherent ejections going 
from the near wall region up to the edge of the boundary layer. These are 


the large structures that recently attractec 
of researchers cited by Pirozzoli et al. [26 


the interest of a large number 
. A more quantitative analysis 


of the effects of the large scale structures has been described by Pirozzoli et 


al. 


261 ] for Couette flows where the turbulent kinetic energy production, in 


the central region, facilitates the formation of the large structures similar to 
those in the boundary layers. 

The transitional region is different for the square bar, the black region in 
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e) f) 

Figure 9: Contours of a), b) u\ and c), d) U 2 at different x 2 — X 3 planes for the cube in 

the top figures and for square bars in the b2®tom at Re = 7500; Aui = 0.05 for positive 
values < 1 colour from blue to red, Aui = 0.005 for values > 1 grey and black for negative 
values; Au 2 = 0.005 yellow negative red positive, in e) and f) contours of W 2 at X 2 = 0.05 
in a Xi — X 3 plane yellow negative and red positive with Au >2 = 1- 




a) b) 


Figure 10: Mean streamwise velocity profile in wall units: a) for the cube ( CU ), b) for the 
square bar (SQ) at three stations (lines) and compared with the (solid symbols) Schlatter 


et al. 


27] data at R T = 1000 and R T = 1400. 


the contour plots in figure [9}o highlights the large separation bubble which 
does not reattach uniformly in agreement with the contours in the second 
plane. The comparison between the contours at the fourth plane in 
figure |9fc and figure EDI stresses large undulations in the outer layer for the two- 
dimensional obstacle, due to the fast growing disturbances promoted by the 
inflectional point in the profiles of u \, located in the black region in figure |9 }d. 
The contours of «2 in the fourth X 2 — x% plane in figure EDI confirm the strong 
interaction between the near wall and the outer regions. This observation is 
corroborated by the wider region where the profiles of u' 2 , for wall bounded 
flows are constant with respect to that for u\ and u' 3 . The contours of o ;2 in 
figure [9f show that the near wall streaks close to reattachment line are short, 
and long in the fully turbulent region. 

To have a quantitative measure of the distance necessary to get the fully 
turbulent regime the profiles in wall units may help. In boundary layers the 
correct averaging is in time and in the spanwise direction. A further averaging 
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in Xi for a dimensionless length equal to 0.1, accounting for approximately 
50 profiles, allows to store less fields. The profiles are calculated at three 
distances from the obstacle: equal to 115 k, 140 k and 165 k. In this range 


1200 > Rg > 1400 therefore the data of Schlatter et al. [27J at R e = 1000 
and Rg = 1400 are considered. Figure [TUh shows that for the cube at the 
first station (red line) the profile of U + does not have a satisfactory log law, 
further downstream the agreement is better, but the profiles fit the log law 
with B = 4.8 instead of B = 5.0. For the square bar a much better agreement 
is observed in Figure flOb where only at X\ = 115 k the correct log law is not 
fitted. The mean velocity profiles approach faster the fully turbulent one, 
however a fully turbulent regime is established when all the statistics have 
the appropriate profile. 

In figure [TT] the profiles of the normal stresses at the same locations of the 
U + in figure 11 reinforce the observation that the fully turbulent condition 
is reached at a shorter distance for the two- than for the three-dimensional 
obstacles. The profiles in the near wall region converge faster implying that 
the near wall structures form closer to the disturbance, and that the large 
structures in the log and in the outer regions are affected by the inlet condi¬ 
tions. The comparison among the three stresses highlights large deviations 
of u' 2 from the correct value, and this fact corroborates our view that, in 
wall bounded flows, this stress accoimts for the changes at the wall, as it 
was shown by Orlandi et al. 23) and for the disturbances affecting the outer 
region. 


Figure 12 shows large differences between the Schlatter et al. (27) results 
and the present ones. The reason is related to the insufficient Li to get 
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Figure 11: Normal stress velocity profiles in wall units: a), c), e) for the cube ( CU ), 
b), d), f) for the square bar ( SQ) at three stations (lines) and compared with the (solid 
symbols) Schlatter et al. 27| data at R T = 1000 and R T = 1400. 













Figure 12: a) Mean streamwise velocity and b), c), d) normal stress velocity profiles in 
wall units for the square bar (SO) at two stations (lines) ( R T « 1600) compared with the 
(solid symbols) Schlatter et al. (27] data at R T = 1000 and R T = 2000. and the Erm & 
Joubert Q (open symbols) experiments at R T = 1000,1500 and R T = 2000. 
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the fully turbulent regime, to the smaller resolution than that in Schlatter 
et al. 27] and to the different type of disturbances. Therefore, a more 
refined simulation, with the parameters previously reported, was performed 
for the square bar, having observed to be the most efficient to promote the 
transition. This simulation is a first attempt to reproduce the most complete 
and documented experiments in turbulent boundary layers by Erm & Joubert 
|7]. The reproduction of the small obstacles used in the experiments requires 
an enormous amount of computational time , therefore in the discussion it 
should be recalled that the height of the tripping device (k/5 = 0.25) is 
greater than that in the experiment (k/8 = 0.1). This difference together 
with the limitations in L 1 affect the wake region. Far from the obstacles 
Rg ~ 1600 is a value intermediate between the Rg = 1000 and Rg = 2000 
by Schlatter et al. [27]. The very good agreement for U + is appreciated 
in figure H2h . Near the wall the stresses (figure fl2b.df are slightly smaller 
than those by Schlatter et al. [2^]; the reason may be related to the different 
resolution. The small oscillations in the outer regions, in large part, are 
due to the large scales produced by the obstacle. In a longer domain these 
oscillations disappear, but to show it a very large CPU is required. These 
oscillations, due to the trip, are not encountered in simulations based on a 
recycling procedure, where these long-lived structures, in the outer region, 
are not continuously introduced. At Rg ~ 700, Erm & Joubert [8| claimed 
that differences with the Spalart 3l) DNS were related to the tripping. In the 
experiments the height was 0.1<5, therefore, with an obstacle of height 0.25J 
it is reasonable to fold in the outer layer differences between the present and 
the results obtained with the recycling procedure. 
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Figure 13: Spectra in kolmogorov units a), c), e) at y + = 25 , b), d), f) at y + = 110 in 
the spanwise direction (lines), compared with isotropic turbulence at R\ = 100 (Jimenez 
et al. Q) open symbols, a) , b) ,c), d) transverse e), f) longitudinal. 
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To demonstrate that the turbulent fields are fully resolved it is worth 
looking at the spectra in Kolmogorov units, if a collapse in the exponential 
range with those of isotropic turbulence occurs it means that a true DNS 
is performed. The agreement in the whole range implies the reproduction 
of the turbulent energy cascade process. The velocity spectra in the span- 
wise direction are compared with the longitudinal and transverse spectra by 


Jimenez et al. 


13] at R\ = 98.7. The spectra at y + = 25 at different lo¬ 


cations, separated by a distance equal to 3.25 starting from x\ = 6.4, are 
presented in figure [131 At the first location, inside the separation bubble, 
the spectra are completely different than the others. In the transitional and 
in the fully turbulent regions a good collapse of the spectra at high /c 3 is ob¬ 
tained, corroborating the previous statement that the small scale adjust in a 
short distance. In agreement with figure 13 the energy containing scales for 
are those adjusting faster. Instead those for w 2 require a longer distance 
to reach the streamwise independence in the fully turbulent regime (indi¬ 
cated by small triangles in figure 13). The spectra for u\ in the log region 
at y + = 110 (the figures on the right of figure ITTil) have a reasonable wide 
inertial range. In addition, a better collapse of the spectra of u 2 is achieved 
(figure mil) . and the tendency towards isotropic turbulence is emphasised in 
figure [T3l l and figure [T3b . From the computational side the last two figures 
are a convincing proof that these simulations are true DNS, having resolved 
scales of size close to the Kolmogorov scales. 
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4. Conclusions 


One common way to promote the transition of laminar boundary layers 
is through tripping devices inserted at a certain distance from the leading 
edge. A large number of experiments were devoted to find a critical Reynolds 
number, but a systematic study of the influence of the shape of the obstacle 
was never attempted, even if it was requested by Klebanoff et al. [l3] - This 
lack of data motivated the present study and from the results achieved it 
was clear that a critical Reynolds number based on the velocity at the top of 
the obstacle (Uk) and on its height k could not be the satisfactory quantity 
marking the transition. In fact it is reasonable to expect that by varying 
the shape of the obstacle, Uk could not drastically change. This has been 
confirmed by the present simulations, but in the mean time it was observed 
that for certain shapes there was transition and for others not. In agreement 
with the previous results in rough channels (Orlandi 20]) it has been observed 
that the transition occurs when u' 2 \k is greater than a threshold value. In the 
channel the periodic conditions lead to an easy determination of the threshold 
value. In boundary layer the determination is more critical for the streamwise 
evolution of u 2 \k- However the numerical simulations demonstrate that, in 
analogy to the case of a normal wall jet where the transition is linked to the 
exit velocity, also for solid obstacles if the amplitude of u' 2 \k is greater than 
a threshold limit the transition occurs. 

The DNS furnishes the vorticity held, and how visualizations of the three 
vorticity components lead to conclude that u >2 is the most important, mod¬ 
ifying the oj 3 sheet generated at the top of the obstacle. The head of the 
horseshoe vortex generated by the obstacle does not contribute to the tran- 
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sition. It has been, also understood that for isolated elements the transition 
is enahnced by geometries which create curved vorticity sheets. To reach a 
fully turbulent conditions the wake generated by isolated obstacles expand 
laterally by the generation of opposite sign vorticity at the wall. The entire 
process can be summarised as the wake instability mentioned by Klebanoff 
et al. [ill, which differs from the inflectional instability occurring in pres¬ 
ence of two-dimensional obstacles. The present DNS demonstrated that the 
latter instability is more efficient to produce a fully turbulent regime, and 
therefore simulations with square bars have been performed at high Re and 
the results have been compared with previous numerical simulations, and 
laboratory experiments. 

The most refined simulation, with a two-dimensional trip, produced mean 
velocity and turbulent stresses profiles in a good agreement with those achieved 
in DNS with recycling procedures or with ad hoc inlet disturbances. The re¬ 
sults suggest that the numerical method here used is an efficient tool to create 
a data base with all the ingredients necessary to understand the differences 
among wall bounded flows in pipe, channels and boundary layers. In the 
simulations ’’hot wires” measuring the three velocity, and vorticity compo¬ 
nents, together with the pressure have been inserted. The temporal date can 
be used to understand the differences on the spectra obtained by the Taylor 
hypothesis and those evaluated in the spanwise directions. The latter spec¬ 
tra were presented as a proof of the reproduction of the physics of turbulent 
flows characterised by an universal exponential range. The temporal signals 
may help to understand the quality of a hot wire, having observed that the 
smallest Kolmogorov scales have been captured in the numerical simulations. 
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The present numerics is not appropriate to reproduce tripping devices of 
smaller height, requiring large computational resources. To achieve this goal, 
with limited resources, modifications of the numerics with non-uniform grids 
in all three directions are necessary. In this way the increased resolution near 
the obstacle allows to create the vorticity sheets driving the transition. For 
incompressible flows this is a difficult task, but it is feasible with non-uniform 
grids in X\ and X2 . Our aim is to introduce these modifications in order to 
reproduce the experiments by Erm & Joubert with a two-dimensional trip 
device with k/5 = 0.1. For three-dimensional obstacles a good clustering can 
be obtained by using the compressible code at low Mach number used by 
Bernardini et al. jo| . In these circumstances the small time steps require 
a much large computational time than that used for the simulations here 
discussed. 
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